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Abstract: We study the influence of the surrounding gas in the dynamics of drop impact on a smooth surface. 
We use an axisymmetric 3D model for which both the gas and the liquid are incompressible; lubrication regime 
applies for the gas film dynamics and the liquid viscosity is neglected. In the absence of surface tension a 
finite time singularity whose properties are analysed is formed and the liquid touches the solid on a circle. 
When surface tension is taken into account, a thin jet emerges from the zone of impact, skating above a 
thin gas layer. The thickness of the air film underneath this jet is always smaller than the mean free path 
in the gas suggesting that the liquid film eventually wets the surface. We finally suggest an aerodynamical 
instability mechanism for the splash. 

PACS numbers: 47.55.D-, 47.55.nd 
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Introduction: Drop impact is present in many surface 
flows and diphasic dynamics. It is crucial to our under- 
standing of atomization, ink-jet printing or deposition 
for instance^, as well as for environmental issues such as 
raindrop impact erosion or aerosol spreading 2 . The gen- 
eral situation involves an almost spherical drop impact- 
ing on a dry or wet solid surface. A splash is observed 
for strong impact conditions, whereas a gentle spread- 
ing of the drop is seen otherwise. On a dry substrate 
as considered further on, the splash is characterized by 
a corolla that detaches from the solid substrate forming 
a corona shape from which droplets can eventually de- 
tach^. This complex dynamics involves various param- 
eters which can highly influence the transition between 
these two regimes: primarily viscous and capillary ef- 
fects are invoked, quantified by the Reynolds and the 
Weber numbers respectively. The impacted surface is 
also important: in particular, the roughness of the sur- 
face can control the splash formation as enhanced with 
textured surfaces^. Recently, the surrounding gas, often 
neglected in these problems, has been shown to be cru- 
cial since the splashing observed at atmospheric pressure 
disappears when the ambient gas pressure is lowered^. 
Despite the few models already proposed to explain this 
striking effect, invoking in particular compressibility of 
the surrounding gas£ or the entrapment of a (compress- 
ible) gas bubble by air cushioning^, a complete under- 
standing for the mechanism of the splash formation is 
still lacking. The aim of this paper is thus to disentan- 
gle the role of the surrounding gas for the drop impact 
on solid surfaces in the limits where the viscosity of the 
liquid of the drop can be neglected and where both fluid 
can be considered incompressible. 

Governing equations and numerical method. We con- 
sider the impact of an incompressible liquid drop of ra- 
dius R, density pi, with vertical velocity V on a solid 
substrate (see figure [J) . The gas is taken incompressible 
of density p g , dynamic viscosity rj with surface tension 




FIG. 1. Sketch of the impacting drop. 



7. Although in the experiments^ the gas is clearly in a 
compressible regime beneath the drop due to the large 
pressures created by the air cushioning, we argue that 
the gas compressibility might not be crucial to under- 
stand the splashing transition since it does not change 
the general structure of the equations. Gravity can be 
neglected since the Froude numbers (Fr — V 2 /gR) are 
always above 10 2 and axisymmetric approximation can 
be safely assumed at short time. We note t = the time 
at which the drop would touch the wall in the absence of 
the surrounding gas. The dominant effect of the gas lies 
in the dynamics of the viscous thin gas film underneath 
the drop. In this situation, the viscous boundary layer in 
the liquid created by the gas shear flow is small and the 
liquid velocity can be approximated by a potential flow 
while the dynamics of the gas layer follows the lubrica- 
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tion equation&ii. Our model is thus similar to the 2D ap- 
proach addressed recently- ^ 11 ' 12 , where a liquid parabola 
impacts a solid surface, but with crucial differences: 1) 
it is 3D-axisymmetric; 2) a spherical drop impacts; 3) we 
use a curvilinear description of the interface which allows 
for the description of the jet. 

Finally, the liquid drop and the gas film dynamics are 
given by the following dimensionless set of equations in 
the cylindrical coordinates (r, z): 
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where tp is the velocity potential in the drop (u(r, z, t) — 
S7(p(r, z,t)), ri is the liquid domain, dn is the boundary of 
the drop and p is the lubrication pressure in the gas. The 
lengths have been rescaled by R, velocities by V, densities 
by pi and the gas pressure by piY 2 . The full interface 
{r(s, t), z(s, t)} is indexed by the curvilinear coordinate 
s and k is the mean curvature of the interface. This set 
of equations introduces the two dimensionless numbers 
of the problem, the Weber and Stokes numbers : 



We 
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and St 



piVR 



Remarkably, general experimental conditions correspond 
to St < 1 as considered further on. The incompress- 
ibility condition Q joint with the Bernoulli equation 
at the drop interface ([I} describes the liquid potential 
flow. This dynamics is coupled with the surrounding gas 
flow through the interface advection equation ([3]) and 
the pressure, given in the thin film by the lubrication 
approximation ([2]). Notice that the lubrication is only 
valid when the air film is thin enough, otherwise a free 
surface flow condition (p = Pq constant pressure) has 
to be applied in the Bernoulli equation. The numerical 
method proceeds as follows : Laplace's equation (j4]) is 
solved using a boundary integral method, the pressure 
is calculated through the lubrication equation ([2])—, the 
interface and the velocity potential are advanced in time 
using the kinematic condition ([3]) and Bernoulli equation 
(HJ) respectively. 

A sequence of snapshots of the drop impact is shown 
on figure[T]for We = 23.7 and St = 1.35 x 1(T 3 . The drop 
deforms as it approaches the wall and a dimple appears 
underneath. Then a quasi-horizontal liquid jet expands 
rapidly. We observe on figure [1] that the liquid never 
touches the wall and the jet "skates" on a thin gas layer! 
This is consistent with the general property of viscous 
film that cannot break-up in a finite-time^. However, 
this is not the case in the absence of surface tension where 
corner like interface can be created as described below. 

Finite time singularity for We = oo: as observed in 
2D & 11 ! 12 , the dynamics exhibits a finite time singularity 
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FIG. 2. Interface profile near the singularity for a 
with zero surface tension and St = 1.35 x 10~ 3 . 
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in the zero surface tension shown on figure 

It corresponds as t — > to to a corner like interface located 
at r c (t) — > ro where the curvature Ko(t) diverges as the 
film thickness ho(t) vanishes. Similarly, the maximum 
pressure po(t) is located in r p (t) ^ r c (t) and diverges 
when t — > to (and r p (t) — > ro), following: 

,-i±0.05 i-2±0.05 

p (xh 2 K cch 

as presented on figure ((3]). The radial position ro(t) 
of ho(t) follows approximately the geometrical intersec- 
tion between the undeformcd drop and the solid wall 
(ro ~ \/2i and ro = 1/V2t). Using the typical gas film 
thickness H* = St 2 ^ 3 by balancing the drop inertia and 
the lubrication pressure^, one can estimate to ~ H* and 
thus r ~ St 1 / 3 (fo ~ St -1 / 3 ) at the singularity, in good 
agreement with the numerical results. 

To understand the properties of this singularity, we 
seek self-similar solutions of the form : h(r,t) = 
h (t)H(R), p(r,t) = p {t)P(R), and <p(r,t) = 
<p {t)$(R,Z), where R = (r - r (t))/l(t) and Z = z/l(t). 
Near the singularity, the spatial variations are seen nu- 
merically slower than (to — t) (see caption of figure 
13]), so that the time derivatives shall be replaced by 
(h f /l(t))d R and equations (11121314ft give at dominant 
order, after dropping the tilde : 
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Taking f = St 1 ^ 3 , two regimes can be distinguished, 
depending on the relevant terms in equation [7J 
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FIG. 3. po and ko as functions of ho- The two dotted lines 

show respectively the h 5 and Hq 2 scalings. The inset shows 
ho(t) and po(t) as function of (to —t) close to the singularity, 
with the fitted scalings ho(t) <x (to — t) 3 ^ 2 and po(t) oc (to — 



FIG. 4. Minimum gas film thickness ho for four different val- 
ues of the Stokes number (St = 3.93x 10~ 4 , 7.86x IO" 4 , 1.57 x 
10~ 3 ,2.36 x IO" 3 ;, as a function of time. The zoom for 
St = 3.93 x IO -4 near the axis shows the evolution for dif- 
ferent Weber number, We = 47.5, 95, 238 and oo from top to 
bottom. 



Regime I : ho -C I: the second term both in equa- 
tion [5] and in the right hand side of equation [7] can be 
neglected, and balancing all the other terms, we obtain: 
ip ~ St~ 1/3 /i , I ~ St~ 2/3 /ig /2 which leads to the ob- 
served numerical scalings: 

Po ~ h and kq ~ St 4//3 /i^~ 2 . 

We thus deduce that this regime corresponds in fact to 
thick gas layer ho » St 4 / 3 . 

Regime II : ho 3> I: here the dominant terms balance 
gives: tp ~ St _5/3 /ig, I ~ St~ 4/3 h 2 and 

Po ~ St~ 2/3 and k ~ St 8/3 /ig 3 . 

Consistently, this regime holds for thin gas layer ho <C 
St 4 / 3 and it has actually never been reached in all the nu- 
merical simulations so far— i!2r— . This analysis suggests 
that the self similar behavior observed in the numerics is 
valid only for large enough h and it predicts that another 
self-similar regime would hold when approaching the sur- 
face closer. In particular, it is not clear whether or not a 
finite time singularity would still exist (as a cusp then). 
Moreover, the lubrication hypothesis would not be valid 
anymore there, since h /l 3> 1, and full Navier-Stokes 
equations should be considered in the gas film. 

Jet formation with surface tension: when adding the 
surface tension, the singularity is regularized since the 
high curvature regions are smoothed by the capillary 
pressure, as illustrated in figure Q where the minimal 
air film thickness is shown for different Weber and Stokes 
numbers. We observe that the dynamics separates from 
the We = oo case when ho becomes small enough while 
the film thickness converges. As the Stokes number de- 
creases, the dimple size decreases, and the jet appears 



earlier with more capillary waves. Moreover, the small 
angle between the tilted jet and the surface varies both 
with the Weber and the Stokes numbers (see figure (J5|)). 
The minimal gas layer thickness can be estimated thanks 
to the two different self-similar regimes exhibited for the 
singular case We = oo. Indeed, balancing the capillary 
term We _1 K with the singular pressure gives for the gas 
layer depending on the regime: 

I : ho ~ St 8 / 9 We- 2 / 3 and II : ho ~ St^We" 1 / 3 

while the two behaviors cross for We ~ St~ 2/3 (with high 
Weber number for regime II). Such dependence is inves- 
tigated on figure [6] where the minimal gas thickness is 
compared with the predicted scalings. We find that the 
two exponents vary following h ~ St ' 9-1 ' We -0 33-0 - 4 
which is in reasonable agreement with regime II. Fi- 
nally, taking the capillary length due to the drop decel- 
eration y/-fRH*/piV 2 ~ St 1/3 We~ 1/2 for the jet thick- 
ness^, mass conservation^ gives for the jet velocity 
Vjet ~ We 1//2 St -2 / 3 , in good qualitative agreement with 
figure ((5|). 

In conclusion, can we explain now the splashing depen- 
dence on the gas pressure observed in the experiments 7 ? 
Since the two control parameters (We and St) do not vary 
when the gas pressure changes, one has in fact to intro- 
duce another physical mechanism. This was somehow the 
main argument to invoke gas compressibility to explain 
the experiments although no clear mechanism is identi- 
fied yet* 7 -'^ 2 ". Here, we propose a purely incompressible 
mechanism which would be still relevant when gas com- 
pressibility will be taken into account? 8 -. Computing the 
minimal gas thickness in the experiments 7 - following the 
scaling laws for h obtained above, we observe that it is 
of the order of a few A, much below the mean free path 
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FIG. 5. Five different computations, with We = 

238; 95; 47.5; 23.7 respectively, from top to bottom and 

St = 1.35 x 1CT 3 . The bottom figure is for We = 23.7 and 
St = 6.29 x 10~ 4 . 



&3 

■■ p 



10" 





~ ~ - - 


■ vF e -i/3 

We- 2 ' 3 ■ 


■ ur 1 






: ' stw '■ 




' io- 2 






io- ;! 






io- 3 





10 2 



117 



FIG. 6. Minimum gas film thickness ho for four different val- 
ues of the Stokes number (St = 3.93 x 10~ 4 , 7.86 x 10~ 4 , 1.57 X 
10~~ 3 ,2.36 x 10~ 3 ), rescaled by St and not-rescaled (inset), as 
a function of the Weber number. 



even at atmospheric pressure (~ 60 nm). Therefore, it 
is reasonable to consider that the liquid film touches the 
solid and the splashing mechanism would then be related 
to the dynamical instability of a rapid thin liquid him 
expanding on a solid surface. Beside the contact angle 



dynamics, we want to emphasize that such rapidly ex- 
panding film is subject to an aerodynamical instability 
similar to the Kelvin-Helmoltz one of a liquid jet in a 
gas environment. If such instability analysis remains to 
be done, it is likely that it will be dependent on the gas 
density, by analogy with the classical results of Squire 1 ^. 

Finally, since the scaling laws for ho are deduced before 
the jet formation, it is important to notice that liquid vis- 
cosity would not drastically affect the main conclusions of 
our calculations although it would change the jet thick- 
ness (then of order St 1,/3 Re _1 / 2 in general much smaller 
than the capillary length calculated before) and the jet 
velocity (Re 1/2 )^^. 
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